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Abstract — The design of modulation schemes for the physical 
layer network-coded three-way wireless relaying scenario is 
considered. The protocol employs two phases: Multiple Access 
(MA) phase and Broadcast (BC) phase with each phase utilizing 
one channel use. For the two-way relaying scenario, it was 
observed by Koike- Akino et al. |4|, that adaptively changing 
the network coding map used at the relay according to the 
channel conditions greatly reduces the impact of multiple access 
interference which occurs at the relay during the MA phase and 
all these network coding maps should satisfy a requirement called 
exclusive law. This paper does the equivalent for the three-way 
relaying scenario. We show that when the three users transmit 
points from the same 4-PSK constellation, every such network 
coding map that satisfies the exclusive law can be represented by 
a Latin Cube of Second Order. The network code map used by 
the relay for the BC phase is explicitly obtained and is aimed at 
reducing the effect of interference at the MA stage. 

I. Background And Preliminaries 

The concept of physical layer network coding has 
attracted a lot of attention in recent times. The idea of 
physical layer network coding for the two-way relay channel 
was first introduced in [1], where the multiple access 
interference occurring at the relay was exploited so that the 
communication between the end nodes can be done using 
a two stage protocol. Information theoretic studies for the 
physical layer network coding scenario were reported in [2], 
[3]. The design principles governing the choice of modulation 
schemes to be used at the nodes for uncoded transmission 
were studied in [4]. An extension for the case when the 
nodes use convolutional codes was done in [5]. A multi-level 
coding scheme for the two-way relaying was proposed in [6]. 

We consider the three-way wireless relaying scenario 
shown in Fig. 1, where three-way data transfer takes place 
among the nodes A, B and C with the help of the relay 
R. It is assumed that the three nodes operate in half-duplex 
mode. The relaying protocol consists of two phases, multiple 
access (MA) phase, consisting of one channel use during 
which A, B and C transmit to R; and broadcast (BC) 
phase, in which R transmits to A, B and C in a single 
channel use. Network Coding is employed at R in such a 
way that A(/B/C) can decode B's and C's(/A's and C's /A's 
and B's) messages, given that A(/B/C) knows its own message. 
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Fig. 1. A three-way relay channel 



For a two-way wireless relay channel, it was observed in 
[4] for 4-PSK, that for uncoded transmission, the network 
coding map used at the relay needs to be changed adaptively 
according to the channel fade coefficient, in order to minimize 
the impact of multiple access interference. In other words, 
the set of all possible channel realizations is quantized into 
a finite number of regions, with a specific network coding 
map giving the best performance in a particular region. It is 
shown in [7] for every M-PSK constellation used by both 
the users, that every such network coding map that satisfies 
the exclusive law is representable as a Latin square and this 
relationship can be used to get the network coding maps 
satisfying the exclusive law. A Latin Square of order M is 
defined to be an M x M array in which each cell contains 
a symbol from %t = {0, 1, ...t— 1} such that each symbol 
occurs at most once in each row and column [8]. 

While most of the research has been done for two-way relay 
channels, some work has been done for the relay channels 
with three or more user nodes as well. Liu and Arapostathis 
[9], proposed a joint network coding and superposition coding 
for three user relay channels, and claim that the results can 
be easily extended to information relaying cases with more 
than three relay nodes. In this scheme, two stage operations 
are required for encoding and decoding, and four channel 
uses are required for the information exchange, three channel 
uses for the MA phase, and one channel use for the BC 
phase, while for our scheme totally two channel uses suffice. 
In [9], the first three channel uses are utilized by each user 
node transmitting its packet to the relay node. Then the relay 



makes two XOR-ed packets and superimposes them together 
for broadcast. The relay makes two XOR-ed packets, the 
packet from the node with the worse channel gain is XOR-ed 
respectively with the other two packets. Pischella and Ruyet 
[10] also discuss the three-way wireless relaying scenario, 
and propose a method of information exchange among the 
users composed of alternate MA and BC phases. The physical 
layer network coding strategy for this relaying protocol is a 
lattice-based coding scheme combined with power control, 
so that the relay receives an integer linear combination of 
the symbols transmitted by the user nodes. This scheme 
also, however, consists of four channel uses. Park and Oh, 
in [11], propose a network coding scheme for the three-way 
relay channels and present a 'Latin square-like condition' for 
the three-way network code. They also discuss schemes in 
order to improve these codes using cell swapping techniques. 
In this work, though Latin Cubes have been suggested as 
being equivalent to the map used by the relay, the number of 
channel uses the scheme uses is five, and the work doesn't 
deal with the channel gains associated with the channels 
explicitly. In [12], authors Jeon et al. adopt an 'opportunistic 
scheduling technique' for physical network coding where 
users in the MA as well as the BC phase are selected on the 
basis of instantaneous SNR using a channel norm criterion 
and a minimum distance criterion and plot graphs to show 
that the proposed scheme outperforms systems without this 
selection. Their approach, however, utilizes six channel uses. 

For our physical layer network coding strategy, we use a 
mathematical structure called a Latin Cube, that has three 
dimensions out of which one is represented along the rows, 
one along the columns, and the third dimension is represented 
along 'files'. In our case, we have A's transmitted symbol 
along the files, B's symbol along the rows, and C's symbol 
along the columns. For our purposes, we define Latin Cubes 
as follows: 

Definition 1: A Latin Cube L of second order of side M 
on the symbols from the set Z t = {0, 1, 2, t — 1} is an 
M x M x M array, in which each cell contains one symbol and 
each symbol occurs at most once in each row, column and file. 

The above definition, is given in [13] with t = M 2 . 

A. Signal Model 

Multiple Access Phase: 

Let S denote the symmetric 4-PSK constellation 
{±1, ± j} as shown in Fig. 2, used at A, B and C. Assume 
that A(/B/C) wants to send a 2-bit binary tuple to B and C(/A 
and C/A and B). Let [i : F 2 , — ► S denote the mapping from bits 
to complex symbols used at A, B and C where F2 = {0, 1}. 
Let xa = ll(sa),xb = h(sb),xc — l^(sc) € S denote 
the complex symbols transmitted by A, B and C respectively, 
where sa,sb,sc G F 2 ,. ^ ^ s assumed that the channel state 
information is not available at the transmitting nodes A, B 
and C during the MA phase. The received signal at R in the 




F 2 XF 2 = {00, 01, 10, 11} 
Z4={0, 1, 2, 3 } 

s= { U-1-J} 



Fig. 2. 4-PSK constellation 



MA phase is given by 

Y R = H A x A + H B x B + H c x c + Z R (1) 

where Ha, Hb and He are the fading coefficients associated 
with the A-R, B-R and C-R link respectively. The additive 
noise Zr is assumed to be CM (0,<7 2 ), where CAf[0,a 2 ) 
denotes the circularly symmetric complex Gaussian random 
variable with variance cr 2 . 

Let Sr (Ha, Hb, He) denote the effective constellations 
seen at the relay during the MA phase channel use, i.e., 

S R (H a ,H b ,H c ) = {H A Xi + H B Xj + HcXklx^xj^Xk E S} . 

Let d m i n (Ha, Hb, He) denote the minimum distance 
between the points in the constellation Sr (Ha, Hb, He) as 
given in @, where S n = S x S x .. x S (ntimes). 
From (fJJ, it is clear that there exists values of 
(H A ,H B ,H C ), for which d min (H A , H B , H c ) = 0. 
Let H = {(HA,HB,H c )eC 3 \d mm (HA,H B ,Hc)=0}. 
The elements of H are called singular fade states. For singular 
fade states, \S R (H A , H B ,H C )\ < 4 3 . 

Definition 2: A fade state (Ha, Hb, He) is defined to be 
a singular fade state for the MA phase of three-way relaying, 
if the cardinality of the signal set Sr(Ha,Hb,Hc) is less 
than 4 3 . Let H denote the set of all singular fade states for 
the three-way data transfer among A, B and C. 

Let (xa,xb,xc) £ S 3 denote the Maximum Likelihood 
(ML) estimate of (xa, xb, xc) at R based on the received 
complex number Yr, i.e., 



(x A ,x B ,x c ) = arg min \\Y R - HX\\ (9) 



where, 

H = [H A H b H c ] 

XA 

X = x B 
. x c 

Broadcast (BC) Phase: 

The received signals at A, B and C during the BC phase 
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d min (H A ,H B ,H c ) = min \H A (x a - x A ) + H B (x B - x B ) + H c (xc — x c ) ] (2) 

( x A , x B, x c),( x ' A , x B< x 'c) eS 

( x A, x B, x c)^( x ' A , x 'b' x 'c'> 

d mVn 3 (Ha,H b ,H c ) = ^to \H A (x A - x' A ) + H B (x B - x' B ) + H C {x c - x' c ) I (3) 

( x ' A , x B> x 'c)£ C 3 

d min (c Ha ' Hb ' H c) = min " " \H A (x A - x' A ) + H B (x B - x' B ) + H c (x C - x' c ) \ . (4) 

(x A ,x B ,x c ),(x' A ,x' B ,x' c )eS , 

M H A' H B. H C(x A ,x B: x c )^M H A. H B. H C (x' a ,x' b ,x'c) ^-^—^ 

M Z1,Z2 (x A ,x B ,xc) 7^ M Z1 ' Z2 (x A , x' B ,x' c ) ,Vx A ,x B ,x' B , xc,x' c 6 S, whenever (x B , xc) ^ (ij, x' c ) (5) 
M Z1 ' Z2 (x A ,x B ,xc) ^ M Z1 ' Z2 (x' A ,x B ,x' c ) yx A ,x' A ,x B ,xc,x' c € S, whenever (x A , xc) ^ (x' A , x' c ) (6) 

M Z1,Z2 (x A ,x B ,x c ) M Z1,Z2 (x' A ,x' B ,x c ) ,Vx A ,x' A ,x B ,x' B ,x c £ S, whenever (x A ,x B ) ^ (x' A ,x' B ) (7) 

(c h ^ h Z' h c,H A ,H B ,H c ) = mh^ \H A (x A - x' A ) + H B (x B - x' B ) + H c {x c - x' c ) \ . (8) 

(x A ,x B ,x c ),(x' A ,x' B ,x' c )£S , 
M h A. h B- h C(x A ,x B ,x c )^M h A. h B< h C(x' A ,x' B ,x' c ) 



'■mm 



are respectively given by, 

Y A = H' A X R + Z Al Y b = H' B X R + Z B , Y c = H' C X R + Zc 

(10) 

where X R = M Ha ' Hb ' Hc ((x a , x b , xc)) G S' is the 
complex number transmitted by R. The fading coefficients 
corresponding to the R-A, R-B and R-C links are given 
by H A , H B , and H c respectively and the additive noises 
Za, Zb and Zc are CN (0, a 2 ) . Depending on the values 
of H A , Hb and He, R chooses a many to one map 
M Ha ' Hb ' Hc : S 3 ->• S' where S' is a signal set of size 
between 4 2 and 4 3 used by R during BC phase. Notice that 
the minimum required size for S is 16, since 4 bits about the 
other two users needs to be conveyed to each of A, B and C. 

The elements in S 3 which are mapped to the same signal 
point in S by the map M h a,h b ,h c ^ t0 f orm a 
cluster. Let {£i, £2, •■, £1} denote the set of all such clusters. 
The formation of clusters is called clustering, denoted by 

qH a ,Hb,Hc^ 

Definition 3: The cluster distance between a pair of clusters 
£i and £j is the minimum among all the distances calculated 
between the points (x A , xb,xc) G £i and (x'a, x'b, x'c) G £j 
in the effective constellation seen at the relay node R, as 
given in d3j above. 

Definition 4: The minimum cluster distance of the 
clustering C Ha ' Hb,Hc is the minimum among all the cluster 
distances, as given in dU at the top of this page. 

In order to ensure that A(/B/C) is able to decode B's and 
C's(/A's and C's /A's and B's) message, the clustering C 
should satisfy the exclusive law, as given in (O, ©, © above. 

The minimum cluster distance determines the performance 
during the MA phase of relaying. The performance during the 
BC phase is determined by the minimum distance of the signal 
set S . For values of (Ha, Hb, He) in the neighborhood 
of the singular fade states, the value of d m i n (C h a,h b ,h c \ 



is greatly reduced, a phenomenon referred to as distance 
shortening [4]. To avoid distance shortening, for each singular 
fade coefficient, a clustering needs to be chosen such that the 
minimum cluster distance at the singular fade state is non 
zero and is also maximized. 

A clustering q h a,h b ,h c j s sa ^ to remove singular fade 
state {H A ,H B ,H C ) G H, if d mm (C^.ife.ifc) > _ For a 
singular fade state {H A ,H B ,H C ) G H, let C^ Ha ^ Hb ' Hc ^ 
denote the clustering which removes the singular fade state 
(Ha, Hb, He) (if there are multiple clusterings which 
remove the same singular fade state (Ha, Hb, He), consider 
a clustering which maximizes the minimum cluster distance). 
Let C H = {cU Ha < Hb > h °» : (Ha,H b ,H c ) G H} denote 
the set of all such clusterings. 

Definition 5: The minimum cluster distance of 
the clustering C hA,hB ' h °, when the fade state 
(Ha, Hb, He) occurs in the MA phase, denoted by 
d m in {C hA ' hB ' hc ,Ha,H b ,H c ), is the minimum among all 
its cluster distances, as given in (0. 

For (Ha,Hb,Hc) £ H, the clustering 
c h a m b ,h c is c hosen to be C {( - hA ' hB ' hc)} , which 
satisfies d mm (C« hA > hB ^\HA,H B ,H c ) > 

dmin (c{( h ' A > h ' B > h 'c)},H A ,H B ,Hc) ,V(h A ,h B ,hc) ^ 
(h' A , h' B , h' c ) 6 T-L. In [7], such clusterings that remove 
singular fade states are obtained with the help of Latin 
Squares while concentrating only on the first minimum 
cluster distance. The clustering used by the relay is indicated 
to A, B and C using overhead bits. 

The contributions of this paper are as follows: 

• Using our proposed scheme, exchange of information in 
the wireless three-way relaying scenario is made possible 
with totally two channel uses. 

• It is shown that if the three users A, B, C transmit points 
from the same 4-PSK constellation, the requirement of 
satisfying the exclusive law is same as the clustering 
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Fig. 3. The mapping observed at the Relay can be viewed as a Latin Cube of Second Order 



being represented by a Latin Cube of second order of 
side 4. To the best of our knowledge, this is the first work 
with only two channel uses for the three-way relaying 
scenario. 

> The singular fade states for the three user case are 
identified. 

« Clusterings that removes these singular fade states are 
obtained, that result in the size of the constellation used 
by the relay node R in the BC phase to lie between 16 
to 23. 

« Simulation results are provided to verify that the adaptive 
clustering as obtained in the paper indeed performs better 
than non-adaptive clustering. 
The remaining content is organized as follows: Section II 
demonstrates how a Latin Cube of Second Order and side 
4 can be utilized to represent the network code for three 
user communication. In Section III the singular fade states 
are specified and in Section IV, clusterings corresponding 
to removal of each singular fade state are obtained using 
Latin Cubes of Second Order. Simulation results are shown 
in Section V. Section VI concludes the paper. 

II. The Exclusive Law and Latin Cubes 

The nodes A, B and C transmit symbols from the same 
constellation, viz., 4-PSK. Our aim is to find the map that 
relay node R should use in order to cluster the 4 3 possibilities 
of (xa,xb,xc) such that the exclusive law (given by (0, 
©, 0) is satisfied. The size of the constellation, or the 
number of clusters of the clustering that the relay utilizes has 
to be at least 16, since each user needs the 4 bit information 
corresponding to the other two users. Consider a 4 x 4 x 4 
array, whose 64 entries are indexed by (xa,xb,xc), i- e - the 
three messages that A, B and C send in the MA phase. Each 
file of this 4x4x4 array, is indexed by a single value of 
xa- Each row (column) of each file is indexed by a value 
of Xb (xc), for a fixed value of xa- Now, a repetition of a 
symbol in a file results in the failure of exclusive law given 
by (0. Consider the 4x4 array with its rows being the 
first(/second/third/forth) rows of the 4x4x4 array. Each 
4x4 array so obtained, corresponds to a single value of xb- 
A repetition of a symbol in this array will result in the failure 
of exclusive law given by (0. Similarly, consider the 4x4 
array with its columns being the first(/second/third/forth) 



columns of the 4x4x4 array. Each 4x4 array so obtained, 
corresponds to a single value of xc- A repetition of a symbol 
in this array will result in the failure of exclusive law given 
by (0. Hence, if the exclusive law needs to be satisfied, 
then the cells of this array should be filled such that the 
4x4x4 array so obtained, is a Latin cube of second order, 
for t > 16 (Definition 1). The clusters are obtained by putting 
together all the tuples (i,j,k),i,j,k e 0,1, ...t — 1 such 
that the entry in the k)-th slot is the same entry from 
Z t . From above, we can say that all the relay clusterings 
that satisfy the mutually exclusive law forms Latin Cubes of 
second order of side 4 for i > 16, when the end nodes use 
4-PSK constellations. Hence, now onwards, we consider the 
network code used by the relay node in the BC phase to be 
a 4 x 4 x 4 array with files(/rows/columns) being indexed 
by the constellation point used by A(/B/C), symbols from 
the set Z4 (Fig. 3). The cells of the array will be filled with 
elements of Z t in such a way, that the resulting array is 
a Latin Cube of Second Order of side 4 and t > 16. Any 
arbitrary but unique symbol from Z t denotes a unique cluster 
of a particular clustering. 

III. SINGULAR FADE SUB SPACES 

We earlier stated in Section II, that a clustering C A' Hb ' G 
is said to remove singular fade state (Ha, Hb, He) £ TL, 
if d min (C Ha < Hb ' Hc ) > 0. Alternatively, removing singular 
fade states for a three-way relay channel can also be defined 
as follows: 

Definition 6: A clustering qH a ,h b m c j s t0 
remove the singular fade state (Ha, Hb, He) 6 H, if 
any two possibilities of the messages sent by the users 

(xa,xb,xc) ,(x' A ,x' B ,x' c ) 6 S 3 that satisfy 

H a xa + H B x B + H c x c = H A x' A + H B x' B + H c x' c 
are placed together in the same cluster by the clustering. 

Definition 7: A set {(xa,xb,xc)} £ S 3 consisting of 
all the possibilities of (xa, xb, xc) that must be placed 
in the same cluster of the clustering used at relay node 
R in the BC phase in order to remove the singular fade 
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state (Ha, Hb, He) is referred to as a Singularity Removal 
Constraint for the fade state (Ha, Hb, He) for three-way 
relaying scenario. 

Let (H A , H B , He) be the fade coefficient in the MA 
phase. The work in [4] and [7] shows that for the two-way 
relaying scenario, the 4 2 possible pairs of symbols from 
4-PSK constellation sent by the two users in the MA 
phase, can be clustered into a clustering dependent on a 
singular fade coefficient, of size 4 or 5 in a manner so 
as to remove this singular fade coefficient. In the case 
of three users, at the end of MA phase, relay receives a 
complex number, given by (Q~|). Instead of R transmitting 
a point from the 4 3 point constellation resulting from all 
the possibilities of (xa, xb, xc), the relay R can choose to 
group these possibilities into clusters represented by a smaller 
constellation. We describe one such clustering in the following. 

Let r denote a singularity removal constraint corre- 
sponding to the singular fade state (Ha, Hb, He) and let 

(x a ,xb,xc),(x'a,x' b ,x' c ) EC. Then, 



xa x'a, xb — x' B and xc — x' c . The singular fade 



H ax a + H B x B + H c x c = H A x' A + H B x' B + H c x' c 
=^Ha(xa - x' A ) + H b (xb - x' B ) + H c (xc - x' c ) = 



X _ 



>(Ha,Hb,H c ) E 



x 



XA 

xb 

xc - x c 

where for a 3 x 1 non-zero vector v over 



(11) 



(v) — {w — (wi,W 2 ,W 3 ) I WiVi + w 2 v 2 + W 3 V 3 = 0} . 

Clearly, (v) is a two-dimensional vector space of C 3 . 
These values of xa , xb , xc , x' A , x' B , x' c E S result in only 
finitely many possibilities for the right-hand side, since S is 
finite. Thus the singular fade states (Ha, Hb, He), which 
are uncountably infinite, are points in a finite number of 
vector subspaces of C 3 . Henceforth, we shall refer to these 
finite number of vector subspaces as the Singular Fade 
Subspaces. More precisely, there are three possibilities of 
singular fade subspaces for the three-way relaying as we 
explain individually in the following three cases. 



Case 

1) 
2) 
3) 
Case 

1) 
2) 
3) 
Case 



1: One of the following subcases arise: 

xa = x' A , xb = x' B and xc ^ x' c 
xa = x' A , xb 7^ x' B and xc = x' c 
xa 7^ x' A , xb = x' B and xc = x' c 
2: One of the following subcases arise: 
xa = x' A , xb x' B and xc ^ x' c 
xa 7^ x' A , xb = x' B and xc ^ x' c 
xa 7^ x' A , xb 7^ x' B and xc = x' c 
3: xa 7^ x' A , xb ^ x' B and xc ^ x' c 



Case 1: Without loss of generality, we discuss 
the third subcase of Case 1, i.e., the case when 



subspace in this case is given by S' = 



xa ~ Xa 





The set of differences of the points of S is given by, 

V = {x t - Xj | Xi, xj e5} = {0, ±1 ± j, ±2j, ±2} . 
Let 2?i = {±1 ± j} and V 2 = {±2j, ±2}. Then, 

D = {0}UDiUP 2 . 

Thus, xa — x' A G V can take eight non-zero values. As 
a result, there are eight total possibilities for the vector 
[xa — x'a, 0, 0]*, where v l denotes the transpose of a vector 
v. Also, each one of ±1 ±j, ±2j and ±2 can be obtained as 
scalar multiples of 1 + j (over C). Thus, 



" 1 + jf " 




" ±l±j - 




" ±2j " 




±2 " 





)■( 





)■( 





)■{ 
























So, 3 only one singular fade subspace for the subcase, viz.. 

J_ 

/ L+J \ 

S' 



1+3 





Similarly, for the other two subcases, there is one singular 
fade subspace, resulting in a total of 3 singular fade subspaces 
for the case. 



The singular fade 



Case 2: Without loss of generality, we discuss 
the third subcase of Case 2, i.e., the case when 
xa ^ x' A , xb ^ x' B and xc = x' c . 

xa — x'a 
xb - x' B 


Here xa — x' A and xb — x' B £ T> can take eight non-zero 
values each. There are therefore, 64 total possibilities for the 
vector [xa — x'a, x B — x' B , 



subspace in this case is given by S" = 



0]*. 



Lemma 1: For the case when xa 



o, X B 



x' B 7^ and xc 



0, for a given vector v 



[xa — x'a, xb — x' b , 0] over T>i UV 2 , there are precisely 
4 or 8 vectors (including v) over T>\ U T> 2 that generate the 
same vector space over C as v. 

Proof: As given in Section II of Q, for the 4-PSK 
constellation S, the difference constellation V = AS = 
{s — s' : s, s' G S} is of the form, 

AS = {0} U |2sm (ttu/4) e° k7T/2 \n oddj 

U {2sin (vm/4) e ^ fc7r / 2+7r / 4 ) |n even} , 
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where 1 < n < 2 and < k < 3. Therefore, we can write, 





XA x' a 




" 2sin^e j < l>1 


V = 


x B - x' B 




2sm^e^ 2 













where fa — kiTi/2 if fcj is odd and = kin/2 + 7r/4 if fcj is 
even. 

A vector w over T>i U X>2 shall generate the same vector 
space over C iff to is a scalar multiple of v, i.e. for some 
complex number re 3 9 G C, 



v = re?® w 



2sm^e^ 2 




re J 



2sin^e^ 3 
2sin^e j ^ 




where for i = 3,4 fa — kiir/2 if ki is odd and 
fa = kiTr/2 + 7r/4 if ki is even. 



Then, 



2sm— e^ 1 = re ]f> x 2sm— e 3 ^ (12) 
4 4 



and 



2sin e JV2 = re J x 2s*n e jq 

4 4 



(13) 



Dividing ( flZb by (fT3l and taking modulus of both sides, we 



get 



(14) 



As shown in Q, this is possible only if k\ = k$ and k 2 = k^. 
Also, from dl~2t and dl~3l we have 



7rfc4 



4 4 

From ( fT~4-b and ( TToT ). we have 

e j{<f>l-<f>2) _ e j(4>3-<f>4.) 



(16) 



Note that here, the LHS is fixed. It therefore suffices to 
compute that for the fixed value of the LHS, the number of 
values that RHS takes. It can be verified, that for a fixed value 
of fa — fa, there are precisely four pair of values of fa and 
fa that result in the same value of fa — fa. We now look at 
the following possibilities: 

Case 1: k\ = k 2 . Then, fc 3 = fc 4 , i.e., there are exactly two 
possibilities for ki and k 2 , viz., k\ = k 2 = 1 and fci = k 2 = 2. 
With two pairs of values for fc 3 and k^ and four pairs of values 
for fa and <p4, we have a total of eight set of values that w 
can take. Hence, in this case, the vector space generated by v 
can be generated by exactly eight other vectors over T>\ \JT>%. 
Case 2: k\ ^ k 2 . Then, ^3 7^ ki, i.e., there is precisely one 
possibility for k\ and k 2 , viz., k\ = k% and k\ = ki. With 
only one possible set of values for fc 3 and k^ and four pairs of 
values for fa and fa, we have a total of four set of values that 



w can take. Hence, in this case, the vector space generated by 
v can be generated by exactly four other vectors over T>iUT> 2 . 

■ 

In this case we end up with 12 singular fade subspaces 
given by the null spaces of the space given on the next page 
in Figure [ill] 

So, 3 12 singular fade subspaces for the subcase each 
being the null space of the above 12 spaces. Similarly, for 
each of the other two subcases, there are 12 singular fade 
subspaces, resulting in a total of 36 singular fade subspaces 
for the case. 

Case 3: In this case, xa 7^ x' A , xb 7^ x' B and xq ^ x' c . The 
singular fade subspace in this case is given by 



(17) 



X A 


- x' A 


X B 


— x' B 


Xc 


~ X C - 



Here each of xa 



' a ■ 



x B — x' B and xq — x' c e T> can 



take eight non-zero values. There are therefore, 512 total 
possibilities for the vector [xa — x' A , xb — x' B , xq — x'c] 1 - 
In this case we end up with 112 singular fade subspaces. We 
now explain this. 

Lemma 2: For the case when xa — x' A ^ 0, xb — 
x' B ^ and xc — x' c ^ 0, for a given vector v = 
[xa — x' A , xb — x' B , xc — x' c ] over T>i U V 2 , there are 
precisely 4 or 8 vectors (including v) over V\ U T> 2 that 
generate the same vector space over C as v. 

Proof: As mentioned in the proof of Lemma 1, from Q, 

V = AS = {0} U |2sm {im/A) e jk7r/2 \n odd j 

U [2sm (tto/4) e ^l 2 ^l^[ n even} , 

where 1 < n < 2 and < k < 3. Therefore, we can write, 

2sm^e^ 1 
2sm^e^ 2 





X A 


— x' A 




V = 


x B 


— x' B 






Xc 


~ x'c - 





2sin^e j <t> 3 



where 
even. 



kiir/2 if ki is odd and fa — kiir/2 + ir/4 if ki is 



A vector w over T>i U T> 2 shall generate the same vector 
space over C iff w is a scalar multiple of v, i.e. for some 
complex number reJ e G C, 



2sin^e j ^ 
2sm^e^ 2 
2sm^e^ 3 



2sm^e^ 5 
2sm^e^ 6 



where for i — 4, 5,6 fa — kiir/2 if ki is odd and 
fa — kiir/2 + 7r/4 if ki is even. 
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Fig. 4. Null Spaces of the Singular Fades Spaces for the case xa ^ x'., xg ^ x'„ and xq = x'„ 



Then, 



2sin^e^ = re' 9 x 2sm^e^\ 
4 4 



Ism e 

4 



>2 _ re je x 2s j n e . 

4 



ttA; 



and, 



2sm — -e 

4 



]<P3 



i6 n • "" fc 6 

re J x isin e 

4 



.7 06 



(18) 
(19) 

(20) 



Dividing (fT8l by (IT~9T > and taking modulus of both sides, we 



get 



sin 



ivk 1 



(21) 

As shown in Q, this is possible only if k\ = k± and k 2 = k§. 
Similarly, dividing ( fT8l by ( 120b and taking modulus of both 
sides, we get 



7rfc 4 



(22) 

ID V I V ^ fD L I Li J " 

As shown in Q, this is possible only if = k± and k 3 = k 6 . 
Also, from ( fT8l and ( fT9] > we have 



sin 



Tiki 
4 g,-( 



sin 



nk 



4 gj'(04- 



5) 



From (f2Tb and d23l) . we have 



Similarly, 



2) = gi( 



gj(0i-03) _ gj(04-0e) 



(23) 



(24) 



(25) 



In d24l and d25l l. the LHS is fixed. It therefore suffices to 
compute the number of values that RHS takes for fixed LHS 
in the two equations. It can be verified, that for a fixed value 
of 4>i, 4>2 and 03, there are precisely four set of values of 
04, 05 and 06 that result in the same value of 0i — 02 and 



0i — 03. We now look at the following possibilities: 
Case 1: k\ = k 2 = k 3 . Then, fct = k$ = fc 6 , i.e., there are 
exactly two possibilities for k^, k$ and ^6, viz., k^ — k$ — 
k§ — 1 and k± — k§ — k$ = 2. With two sets of values for 
&4, fcs and fcg and four sets of values for 04, 05 and 06, we 
have a total of eight set of values that w can take. Hence, in 
this case, the vector space generated by v can be generated by 
exactly eight other vectors over T>\ UV 2 - Case 2: Atleast one 
of fci 7^ k 2 , k\ ^ k$, or k 2 ^ k%. Then, k^ k$, k^ ^ /c6, 
or /c5 ^ ke, so that, using ( EH and (1221 1. there is precisely 
one possibility for k^, k$, ke, viz., k^ — ki, k$ = k 2 and 
k§ — k$. With only one possible set of values for k^, k$ and 
k§ and four sets of values for 04, 05 and 06, we have a total 
of four set of values that w can take. Hence, in this case, the 
vector space generated by v can be generated by exactly four 
other vectors over T>\ UT> 2 . 

■ 

Since all of xa ~ x' A , xb — x' B and xc — x' c £ T> are 
non-zero, we can say that 

xa — x'ai x b — x'b an( l x c — x'c £ T>i U T> 2 . 
As a result, we have the following three subcases: 

1) One of xa — x' A , xb — x' B and xc — x' c £ T>i 

2) Two of xa — x' A , xb — x' B and xc — x' c £ Vi 

3) All of xa — x' A , xb — x' B and xc — x' c £ 2?i 
We deal with each one of the subcases one-by-one. 

Subcase 1: One of x a — x' A , xb — x' B and xc — x' c £ T>\. 
Without loss of generality, we assume that xa — x' A £ T>\ 



and xb — x' B , xc 



£ V 2 . The singular fade subspace 



for the case is given by ( [17) . There are 64 possibilities for 
the vector v' = [x A — x' A , x B — x' B , xc — x' c ] . But some 
of the possibilities may generate the same vector space over 
C. There are precisely 4 vectors of length 3 over T>\ U T> 2 , 
the {±1, ±j} scalar multiples of the vector. As a result, 
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the case xa — x' A G T>\ and xb — x' B , xc ~ x' c 6 T>2 
leads to 16 singular fade subspaces as shown in Figure [HI] 
The same holds for the case when xb — x' B G T>\ and 
xa — x'a, xc — x' c G T>2, or when xc — x' c G T>\ and 
— x' A , xb — x' B G 2?2- This subcase therefore results in 
48 singular fade subspaces. 

Subcase 2: Two of xa — x' A , xb — x' B and xc — x' c G T>\. 
Without loss of generality, we assume that xa — x' a ,x b —x' b G 
T>i and xc — x' c G X>2- The singular fade subspace for the 
case is given by (TTTb and again there are 64 possibilities 
for the vector v" — [xa — x' A , x B — x' B , xc — x'cf ■ F° r 
a given v", possibilities of other 3 length vectors over 
T>i U T>2 that generate the same vector space over C are 
the {±l,±j} scalar multiples of v". As a result, the case 
xa — x' A G T>i and x B — x' B , xc — x' c G T>2 also leads to 
16 singular fade subspaces as shown for this case in Figure 
(Ull). The same holds for the case when x B — x' B G T>2 and 



Xc 



G T>\, or when xa 



G T>2 and 



x B - x B , x c - 
fade subspaces. 



G T>i resulting in therefore 48 singular 



x'a, x b — x' B and xc ~ x[ 



c e 

vector 



Subcase 3: All of — x^, 
2?i. There are 64 possibilities for the 
[xa — x' A , x B — x' B , xc — x'fjf over T>\. For a given 
such vector, possibilities of other 3 length vectors over 
D\ U T>2 that generate the same vector space over C are the 
{±l,±j, ±l±j} scalar multiples of the vector. The case 
xa — x' A , x B — x' B , xc — x'c G T>\ leads to a total of 16 
singular fade subspaces as shown in Figure (11 lib . 

The three cases result in a total of 3+36+48+48+16=151 
singular fade subspaces. We now discuss how these singular 
fade subspaces can be removed using Latin Cubes of Second 
order. 

IV. Removing singular fade subspaces and 
Constrained Latin Cubes 

In the previous section, we classify the set of singular fade 
subspaces into three cases. We now cluster the possibilities 
of (xa,x b ,xc) into a clustering using Latin Cubes. This 
clustering is represented by a constellation given by S', which 
is utilized by the relay node R in the BC phase. The objective 
is to minimize the size of this constellation used by R. 

The clustering to be used at R, first constrains the 
possibilities of (xa, x b , xc) received at the MA phase, with 
the objective of removing the singular fade subspaces, and 
fills the entries of a 4 x 4 x 4 array representing the map 
to be used at the relay using these constraints, and then 
completes the partially empty array obtained to form a Latin 
cube of second order. The mapping to be used at R can be 
obtained from the complete Latin cube keeping in mind the 
equivalence between the relay map with the Latin Cube of 
second order as shown in Section III. 



X A 


— x' A 


xb 


— x' B 


xc 


~ X C - 



In order to to obtain the constraints on the 4x4x4 
array representing the map at the relay node R during BC 
phase for a singular fade state, we utilize the vectors of 
differences, viz., [xa — x' a x b — x' b , xc — x' c ] contributing 
to that particular singular fade state. During MA phase for 
the three-way relaying scenario, nodes A, B and C transmit 
to the relay R. As shown in the previous section, there is 
a total of 151 singular fade subspaces. Let (fiA, h B , he) 
denote a point in one of the 151 singular fade subspaces. 
Then, there exist (xa, xb,xc) , (x'a, x' b , x' c ) G S 3 that satisfy 
Haxa + h B x B + h c x c — fiAx' A + h B x' B + h c x' c . In 
order to remove the singular fade state (Ha, h B , he), the pair 
(xa, x b ,xc) , (x' A , x' B , x' c ) must be kept in the same cluster 
in the clustering. For instance, if 



(h A ,h B ,h c ) G 



then the pair (xa, x b , xc) , (x' A , x' B ,x' c ) must be kept in the 
same cluster in the clustering, i.e., the entry corresponding 
to (xa,x b ,xc) in the 4x4x4 array must be the same 
as the entry corresponding to (x' A ,x' B ,x' c ). Similarly, every 
other such pair in S 3 contributing to the same singular fade 
subspace must be kept in the same cluster. Apart from all such 
pairs in S 3 being kept in the same cluster of the clustering, 
in order to remove this particular fade state, there are no 
other constraints. Upon filling up the 4x4x4 array with 
these constraints, the entire Latin Cube does not get filled up, 
but can be completed as we show later in this section. It is 
important to note that, this clustering cannot be utilized to 
remove the singular fade subspaces of Case 1 and Case 2 of 
the previous section, as shown in the following lemma. 

Lemma 3: The clustering map used at the relay node R 
cannot remove a singular fade state corresponding to the Case 
1 and Case 2 of the previous section and simultaneously satisfy 
the mutually exclusive law. 

xa x a 



Proof: Let S' 



xb 



1 B 



be a singular 



x c -x c 

fade state for Case 1, and without loss of generality, assume 
that xb = x' B and xc — x' c . Then, in order to remove S', 
(xa,x b ,xc) and (x' A , xb, xc) must be kept in the same 
cluster. This is impossible, since this clearly violates the 
mutually exclusive law as if the pair is placed in the same 
cluster, the users B and C will not be able to distinguish 
between the messages xa and x' A sent by the user A. 

Let S' be a singular fade state for Case 2, and without loss 
of generality, assume that xc = x' c . Then, in order to remove 
S', (xa,xb,xc) and (x' A , x' B , xc) must be kept in the same 
cluster. This also clearly violates the mutually exclusive law 
since if the pair is placed in the same cluster, the user C will 
not be able to distinguish between the messages xa, x b and 
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Fig. 5. Null Spaces of the Singular Fades Spaces for the case xa — x' A € T>i and x B — x' B , xc — x' c € T> 2. 
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Fig. 6. Null Spaces of the Singular Fades Spaces for the case xa — x' A , x B — x' B G V\ and xq — x' c £T>2. 



messages x' A , x' B sent by the users A and B. ■ 
The singular fade subspaces given in Case 1 ans Case 2 of the 
previous section, whose harmful effects cannot be removed by 
a proper choice of the clustering are referred to as the non- 
removable singular fade subspaces. The rest of the singular 
fade subspaces, given in Case 3 of the previous section, are 
referred as the removable singular fade subspaces. 

We now illustrate the removal of a Case 3 singular fade 
state with the help of examples. 

Example 1: Let the singular fade subspace be 
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Consider the first vector [1 + j, 2j, — 2j]*. Here, 1 + j 



can be obtained as a difference of xa = 1 and x' A = —j or as 
a difference of xa = j and x' A = —1; 2j can be obtained as a 
difference of xb — j and x' B = — j; and — 2j can be obtained 
as a difference of xc = —j and x' c — j. Thus, the entries 
corresponding to (1, j, — j) and (— j, must be the same 

and the entries corresponding to (j, j, — j) and (— 1, — 
must be the same in the 4x4x4 array representing the 
clustering, i.e., entries (0, 1, 3) and (3, 3, 1) must be the same 
and entries (1, 1,3) and (2, 3, 1) must be the same. 

The second vector [—1 — j, — 2j, 2jf where, —1 — j can 
be obtained as a difference of xa = —1 and x' A = j or as a 
difference of xa = —j and x' A = 1; — 2j can be obtained as 
a difference of xb = — j and x' B = j; and 2j can be obtained 
as a difference of xc — j and x' c = — j. Thus, the entries 
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corresponding to (—1, and (j, j, —j) must be the same 

and the entries corresponding to and 
must be the same in the 4x4x4 array representing the 
clustering, i.e., entries (2, 3, 1) and (1, 1, 3) must be the same 
and entries (3, 3, 1) and (0, 1, 3) must be the same. Similarly, 
the constraints resulting from the third and forth vector ca be 
obtained. The constrained array is shown in Figure 8. 

We choose to fill the constrained array to form a Latin cube 
of second order as shown in adjoining Algorithm 1, which 
simply fills the empty cells of the partially filled 4x4x4 
array with £j, i > 1 in increasing order of i keeping in mind 
that the resulting array must be a Latin cube of second order. 
The top-most and the left-most empty cell in the earliest file 
is filled at every iteration. The completed Latin cube is shown 
in Figure 9. 

Example 2: Consider another example for which the singu- 



lar fade subspace is given by 
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first vector is 




- 1 




A 


Here, 


1 + 



can be obtained as a difference of xa = 1 and x' A = —j 
or as a difference of xa = j and x' A = —1; — 1 — j 
can be obtained as a difference of xc = —1 and 
x' c = j or as a difference of xq = —j and x' c = 1. 
Thus, the entries corresponding to {(1, 1, — 1), (— j, — j, j)}, 
{(1, l, -j), (-j, -j, 1)},{(1, j, -1), (-j, -l, j)},{(l, j, -j), (-j, -l, 1)}, 
{(j, l, -1), (-1, -j, 1, -i), (-1, -j, J. -l). (-1. 

-j), (-1, -l, 1)} must be the same in the 4 x 4 x 4 array 
representing the clustering, i.e., entries {(0, 0, 2), (3, 3, 1)}, 
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Fig. 8. Constraints for Example 1, with B's transmitted symbols along the rows and C's transmitted symbols along the columns 
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Fig. 9. Latin Cube for Example 1, with B's transmitted symbols along the rows and C's transmitted symbols along the columns 



Algorithm 1: Obtaining the Latin Cube of Second Order 
from the constrained 4x4x4 array 

Input: The constrained 4x4x4 array 

Output: A Latin Cube of Second Order representing the 
clustering map at the relay 

1 Start with the constrained 4x4x4 array 

2 Initialize all empty cells of A" to 

3 Let y denote the 4 x 16 matrix obtained by 
concatenating X row-wise, and let Z denote the 16 x 4 
matrix obtained by concatenating X column-wise 



4 The (i,j,k) th element of X is the i th file, the j th row 
and the k th column cell, 
for 1 < i < 4 do 
for 1 < j < 4 do 
for 1 < k < 4 do 

if cell [i, j, k) of X is NULL then 
Initialize c=l 

if C c does not occur in the i th file of X, 
the j th row of y and the k th column of 
Z then 

replace at cell k) of X with 

C c ; 

replace y with the 4 x 16 matrix 
obtained by concatenating X 
row-wise, and Z by the 16 x 4 matrix 
obtained by concatenating X 
column-wise; 

else 

| c=c+l; 
end 



it 
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13 
14 
15 
16 
17 
18 

19 end 



end 



end 



end 



{(0,0, 3), (3, 3,0)}, {(0,1, 2), (3, 2,1)}, {(0, 1, 3), (3, 2, 0)}, 
{(1,0, 2), (2, 3,1)}, {(1,0, 3), (2, 3,0)}, {(1, 1, 2), (2, 2, 1)}, 
{(1, 1, 3), (2, 2, 0)} must be the same. Similarly the other 
constraints can be obtained. The constrained array is shown 
in Figure [To] and the clustering is as shown in Figure [TT] 

For each of the 112 possibilities of singular fade subspaces 



of Case 3 of the previous section, a clustering of size between 
16 to 23 can be achieved by first constraining the array 
representing the relay map in order to remove the singular 
fade state and then completing the constrained array using 
the provided algorithm. 

V. SIMULATION RESULTS 

The proposed scheme is based on the removal of singular 
fade states for the three-way relaying sccenario. A minimum 
cluster distance greater than zero is ensured for all the fade 
states, excluding a subset of singular fade states, which are 
shown to be non-removable. It is attempted to ensure that in 
the given scenario, the number of clusters in the clustering, 
which is the same as the size of the signal set used during 
the BC phase is minimized. Simulation results presented in 
this section identify some cases where the proposed scheme 
outperforms the naive approach that uses the same map for 
all fade states and vice verse All the simulation results shown 
in this section are for the case when the end nodes use 4-PSK 
signal set. The simulation results for the end to end BER as 
a function of SNR are presented in this section for different 
fading scenarios. 

Consider the case when Ha, Hb, He, H' A , H' B and H' c 
are distributed according to Rician distribution and channel 
variances equal to dB. The SNR vs BER curve for this 
case, for a frame length of 256 bits. The plots for the cases 
with a Rician Factors of 5 dB, 10 dB, 15 dB and 20 dB are 
as shown in Fig. 13, Fig. 14, Fig. 15 and Fig. 16 respectively. 
The figures show the SNR vs bit-error-rate curves for the 
following schemes: adaptive clustering given in the paper, 
and non-adaptive clustering. For non-adaptive clustering, the 
relay node uses the same map given by Figure [12] for all the 
channel realisations, we refer to this map as the non-adaptive 
map. It can be seen from Fig. [14] that the schemes based 
on the adaptive clustering relaying perform better than the 
schemes based on non-adaptive clustering at low SNR, since 
adaptive clustering removes 112 singular fade states. 

It can be seen from the simulation results presented that 
the non-adaptive network coding performs better than the 
proposed scheme above a certain SNR when there is a 
dominant line of sight component, as in the case of Rician 
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Fig. 10. Constraints for Example 1, with B's transmitted symbols along the rows and C's transmitted symbols along the columns 
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Fig. 1 1 . Clustering for Example 2 




fading scenario. The reason for this is as follows: The end 
to end SNR vs BER as well as the throughput performance 
depend on the performance during the MA phase as well 
as the BC phase. As the line of sight component becomes 
more and more predominant, the performance during the BC 
phase gets better and better, but the effect of multiple access 
interference which occurs during the MA phase remains the 
same. Hence, for the cases when line of sight component is 
predominant, the performance degradation due to the MA 
interference predominates over the degradation occurring 
during the BC phase. The case of non-adaptive network 
coding, the relay utilizes a constellation of least possible 
size, whereas in adaptive network coding, the relay attempts 
at removing singular fade states, thereby optimizing the 
performance during MA phase to the fullest extent possible, 
at the cost of degraded performance during the BC phase. 
Hence, the non-adaptive network coding performs better than 
the proposed scheme at high SNR. 




Fig. 15. SNR vs ber curves for different schemes for 4-PSK signal set when 
the Rician Factors is 15 dB 
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Fig. 12. Non-Adaptive map 




Fig. 16. SNR vs ber curves for different schemes for 4-PSK signal set when 
the Rician Factors is 20 dB 

VI. Conclusion 

Our paper deals with the three-way wireless relaying 
scenario, assuming that the three nodes operate in half-duplex 
mode and that they transmit points from the same 4-PSK 
constellation. It is shown that it is possible for information 
exchange to take place using just two channels uses, unlike 
the other work done for the case, to the best of our knowledge. 
The Relay node clusters the 4 3 possible transmitted tuples 
(xa, xb, %c) into various clusters such that the exclusive 
law is satisfied. This necessary requirement of satisfying 
the exclusive law is shown to be the same as the clustering 
being represented by a Latin Cube of second order. Using the 
proposed schemes, not only is the exchange of information 
between the three nodes made possible using three channel 
uses, the size of the resulting constellation used by the relay 
node R in the BC phase is reduced from 4 3 to lie between 16 
to 23. Note that we do not claim that the size of the clustering 
utilizing modified clustering is the best that can be achieved, 
since our method of filling the Latin Cube of Second Order of 
side 4 may not be the most optimal process of doing so, and 
it might be possible to fill the array with less than 23 symbols. 
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